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Abstract 

As a toy model of a gapped system, we investigate the entanglement entropy of a 
massive scalar field in 1+1 dimensions at nonzero temperature. In a small mass m and 
temperature T limit, we put upper and lower bounds on the two largest eigenvalues 
of the covariance matrix used to compute the entanglement entropy. We argue that 
the entanglement entropy has e~ m l T scaling in the limit m <C T. We comment on 
the relation between our work and the Ryu-Takayanagi proposal for computing the 
entanglement entropy holographically. 



1 Introduction 



The notion of entanglement entropy (and more generally quantum entanglement) looms 
large in theoretical physics today. Entanglement entropy may be a good order parameter 
for topological phase transitions in condensed matter systems. For conformal field theories 
in 1+1 dimensions, numerical computation of the entanglement entropy provides a rapid 
way to calculate the central charge c. In relativistic field theories more generally, certain 
special kinds of entanglement entropy show monotonicity properties under renormalization 
group flow [H12]- See [3j and [I] for reviews. 

To compute the entanglement entropy for a quantum mechanical system, we must first 
divide the associated Hilbert space up into two pieces. Usually, the division is made with 
respect to spatial regions A and complement A = B. We find the reduced density matrix 
PA = tr# p by tracing over the degrees of freedom in B. Finally, the entanglement entropy 
is defined to be 

S = — tr pa log pa ■ (1) 

It is surprising that even for what many consider to be the simplest field theoretic system 
- a massive scalar field in 1+1 dimensions - the entanglement entropy has thus far been 
computed analytically only in certain limits. In the limit m = 0, one can use results from 
conformal field theory 016]. In particular, for the massless scalar field on the cylinder 
KxS 1 where R is interpreted as the time direction, one has 

1 / L . tt£\ 
5= 3 log^sm T j+c , (2) 

where L is the circumference of the S 1 , £ is the length of the interval, e is a UV regulator 
and Co is a constant that depends on the regulation scheme. (In fact, for the massless scalar, 
there is an additional IR divergence, and cq depends also on an IR cutoff.) Reinterpreting 
S 1 as a Euclidean time direction, one obtains a result at nonzero temperature T = 1//3 for 
the scalar on M. 

When m + 1 for the scalar field on M 2 , Huerta and Casini [7] have shown that the entan- 
glement entropy can be computed from the solution to a certain Painleve equation. Their 
work allows analytic access to the small and large mass limits. For mi <C 1, one obtains 

_ 1, I 1, Aog(me)\ 

S ~ - log - + - log * , 4 

3 e 2 \\og(mi) J 
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while for mi S> 1, one finds instead exponential suppression 

S ~±-pL e -^. (5) 
16 V mi v ' 

Ideally, one would like to understand the case where m, T, and 1/L are all nonzero. 
Numerically, the entanglement entropy can be computed with ease using a generalization [8j 
of the procedure introduced by Srednicki [9j. One realizes the scalar field as the continuum 
limit of an iV-site harmonic chain. For such a chain, one introduces two point functions 
{4>i4>j) and i^iTTj) of the oscillator positions and conjugate momenta respectively. Restricting 
now to an interval ne = i < L where 1 < k < n, one constructs the n x n matrix 

n 

C 2 = J2(<t>i^k)(^j) ■ (6) 

k=l 

The entanglement entropy is then 

S = tr [(C + 1/2) ln(C + 1/2) - (C - 1/2) ln(C - 1/2)] . (7) 

To our knowledge, this quantity has not been computed analytically for the real scalar field 
with two or more of the quantities m, T, and 1/L nonzero. Happily, with today's desktop 
computers, it is relatively quick to diagonalize C numerically for N ~ 10 3 . Ref. jlO] provides 
a numerical analysis of the harmonic chain using this approach. 

In this paper, we take some steps toward an analytic understanding of the eigenvalues of 
C. As noted in |10j . the parity operator P commutes with C where parity here is a reflection 
of the circle S 1 with respect to the midpoint of the interval. Thus, one may divide C into 
even and odd parity blocks C e and C a . We compute the two partial traces trC 2 and trC 2 
in the limit m, T -C 1/L. As the spectrum of C 2 is bounded below by 1/4, these traces give 
us upper bounds on the two largest eigenvalues of C. A variational approach gives a lower 
bound to the largest (parity even) eigenvalue. These bounds in turn give us some intuition 
for the m, T, and L dependence of the entanglement entropy in the limit m,T <C 1/L. 

The original motivation for this project came from our interest in the Ryu-Takayanagi 
proposal [IT] for computing the entanglement entropy of field theories with dual holographic 
classical gravity descriptions. Given two complementary regions A and B in the field theory, 
the Ryu-Takayanagi proposal associates a nonzero Sa~Sb to gravity descriptions with black 
holes, while in the absence of such defects Sa = Sb- In the dual field theory, the existence 
of a black hole typically implies deconfined gauge theory degrees of freedom [121 03] • We 
may contrast this result with the quantum mechanical point of view where at T = 0, the 
density matrix is constructed from a pure state. (We are assuming the existence of a unique 
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ground state.) It follows from a Schmidt decomposition of the Hilbert space that for pure 
states Sa = Sb (see for example [3]). However, at any nonzero temperature, regardless of 
the presence of deconfined degrees of freedom, the density matrix is not constructed from a 
pure state and one would generically expect Sa 7^ Sb- As gauge theories are more difficult 
to study than the free scalar field and as the entanglement entropy of the free scalar field 
has not yet been completely understood, our toy model of confinement in this paper is a 
1+1 dimensional massive scalar field on a circle at T > 0. Morally, the regime T <C m can 
be thought of as "confining" . One of our results is that in this regime, the entanglement 
entropy difference does not vanish but rather scales aa]] 

S A - S B ~ e~ mtT - 



2 From the Harmonic Chain to the Scalar Field 

Consider the Hamiltonian for a real free massive scalar field on a circle of circumference L 
at T > 0: 

H=] 2 jdx [tt(x) 2 + {d x 4>{x)f + m 2 0(x) 2 ] . (8) 
We discretize the circle into N points where L = Ne: 

I N 

where vr(je) = Wj/e but <j>{je) = <j>j. The thermal density matrix can be written in terms of 
H in the standard way: 

P -H/T 

(10) 



tr(e-^/ T ) 

and expectation values are defined via (X) = tr(pX). A short calculation yields the two 
point functions of the oscillator positions <j>j and their conjugate momenta 7Tj-: 

f N-l 

'3W - 



i-V^-coth(^)cosf 27r(j -^ , (11) 

a=0 a v 7 

fa**) = o^r ^ euJa coth iw) cos k ) ' ( 12) 

n— n V / 



a=0 

where 



2 2 , 4 ■ 2 71-0 
uj a = m +^sm — 



1 After finishing this work, we became aware of ref. |14) where the same exponential behavior was found 
for a "renormalized thermal entropy" similar in some respects to the entanglement entropy we study here. 
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From eqs. ^ and ([7]), we may compute the entanglement entropy from the matrix C 2 = 
(tttt) ■ where the two point functions are now restricted to the interval A: —s < j,k < s. 
In terms of n, we have 2s + 1 = n. 

The Hamiltonian H is a set of ./V coupled harmonic oscillators. Diagonalizing the Hamil- 
tonian, one finds H = ^2 a u) a bab a where [b a , b' b ] = 5 a b- Surprisingly for a free scalar field, 
the reduced density matrix pa ~ e~ HA can be written in terms of a similar entanglement 
Hamiltonian Ha = X^fe e kb\bk ( see f° r example jl]). Moreover, there is a one-to-one corre- 
spondence between eigenvalues A^ of C 2 and the energies e^: 

X k = I coth 2 | . (13) 

As the 6^ are real, we conclude that > 1/4. 



3 Taking Traces 

For a region —s<k<s, the matrix C 2 commutes with the parity operator which sends 
k — ^ —k. Thus, we can decompose C 2 into even and odd parity pieces, C 2 = C 2 + C 2 . The 
numerics suggest that for small masses (mL < 1) and low temperatures (TL <C 1), the 
matrix C 2 has only a handful of eigenvalues which are significantly different from 1 /4. The 
larger of these eigenvalues corresponds to an eigenvector with even parity, while the second 
largest has odd parity. We approximate these eigenvalues by computing tvC 2 and trC 2 . 
We find in the even sector that 

' 4N sin(-7rr) s 



trC" 



1 



m 



coth ( 
2vrmL \2T 



7 + In 



7T 



r ,/m 
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+ - 
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AN sin(7rr) 



coth(^) [Li 3 (e 27r -) + Li 3 ( e 



7T 7T 

2izir\ i t ; / — 2mr- 



32vr 3 """" \2TJ 
+0{(mL) 2 , e- 2lT / TL ,log N/N) 



7 + In 
2C(3)] 



AN sin(-7rr) 
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(14) 



and that in the odd sector 



trCf 



s + 



12 



3 1/ , 47Vsin(7rr) 

2^ + 2^V- 1 + ln —^ 



+0((mL) 2 ,e- 27T/TL ,logN/N) , (15) 
where r = i/L. We make some brief remarks about how these traces were computed below. 
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Because of the relation A& = \ coth 2 (efc/2) between the entanglement spectrum and the 
eigenvalues of C 2 , we know that the eigenvalues of C 2 are bounded below by 1/4. The 
largest even eigenvalue A e and odd eigenvalue A G are thus bounded above by 



K < tie: 



Xo < trC 2 



s 

4 ' 
s- 1 



(16) 
(17) 



We can also put a lower bound on A e by using the variational principle and a "trial wave 
function". In this case, we use a constant trial wave function, ip e = (1,1,..., l)/y/n. The 
expectation value then provides a lower bound: 



A e > mC 2 M e ) 
1 



(18) 



/ m \ 



coth , 
2-irmL \2T J 



8vr 3 r 



7 + In 



7 + In 
AN sin(7rr) 



4A^sin(-7rr) 



7T 



+ 



12 



7T 



[Li 2 (e 2 ^)-Li 2 (e- 



■2mr\ 



T 



+ 0(mL,e- 2n/TL , log N/N) . 



Figure [T] demonstrates that our upper and lower bounds provide relatively good estimates 
of the two largest eigenvalues at T = 0. We could try to produce an analytic lower bound 
on A D by similar methods. However, simple trial wave functions such as {ip )j ~ sin(-7rj'/A^) 
or (ip )j ~ J do not seem to give strong lower bounds numerically and are harder to work 
with analytically than the constant trial wave function used above in the even case. 

The zero mode a = terms in (4>4>) and (tttt) have a large influence on the structure 
of these traces in our m, T 1/L limit. As these zero modes have even parity, they 
do not contribute to C 2 . For example, note that trC 2 = 0(l/mL) is much larger than 
trC 2 = O(l) because the zero mode a = term in (</></>) is 0(l/mL) but only contributes 
to the even sector of C 2 . Also note that only trC 2 depends on T. The reason is that 
coth(w a /2T) ?s 1 up to exponentially suppressed terms except when a = 0. 

Another interesting feature of these traces is their behavior under the exchange of the 
interval A with its complement B. By translation invariance, this exchange can be im- 
plemented by sending r — > 1 — r. At T = 0, both trC 2 and trCf are invariant under 
this transformation. This invariance is expected in order to guarantee that Sa = Sb- For 
the breaking of this symmetry is due entirely to the r 2 csch 2 (m/2T) term in tr C 2 . 
This symmetry breaking term comes from multiplying the a = zero modes in (4>(p) and 
{tttt} together. Figure [2] demonstrates that trCf gives a remarkably good estimate of the 
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Figure 1: The largest (a) and second largest (b) eigenvalue of C 2 plotted against the interval 
length for mL = 1/10, T = 0, and N = 1000. The points are numerically computed. The 



curves above the points are the analytic upper bounds (16) and (17) computed from the 



traces. The solid curve below the points on the left is the lower bound (18) computed from 
the variational principle. 



temperature dependence of the largest eigenvalues for regions A and B, and also for their 
difference. 



We should say a few words about the lengthy computation performed to obtain (14), 
(15), and (18). Consider first the 0(1 /mL) contribution to tvC 2 : 



tr C " = 2mT C ° th (If) N) + ° {mL) ° + °^' TL ) > 



where 



N 



_ 1 y> sin 2 (7rtm/iV) 
Jin,ivj _ N 2_^ sin(7ra/A n N 

a=l x ' ' 



n 



(19) 



(20) 



We want to evaluate this sum in the continuum limit where n and N are both large but 
r = n/N is held fixed between zero and one. Replacing the sum over j by an integral 
introduces unacceptably large errors because of the divergence at j = 1/2. Instead, we 
compute a related integral that does not have this divergence: 

1 



f(n,N) 



cot7r(x - 1/2A^) , ^ 

i/n\ V ' ' w(x-l/2N) 

n 



dx 



N 



1 

7T 



^ n(j - 1/2) 
In ( v — | + 7 



7T 



+ 0(1/N 2 



(21) 
(22) 



Calculating the 0(mL)° and 0(mL) terms is a more complicated enterprise. As men- 
tioned already above, one contribution to C? comes from multiplying the zero modes in 



6 




Figure 2: The largest eigenvalue of C 2 as a function of temperature for mL = 1/50: a) 
i/L = 1/5; b) i/L = 4/5; c) the difference between the two for a lattice with N = 200. The 
points are numerical while the curve is the upper bound computed from tr C 2 . 
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and (tttt) together and yields r 2 csch 2 (m/2T). The remaining order one pieces can be 
computed from the matrix C 2 with the zero modes removed in the limit m = = T: 

~ 2 i v^ 1 sin ^ 2nla 27rlb 27T i a 27rkb /« Q \ 
(C p )jk = —tt^ > > cos — cos —— cos Ar cos — . (23) 

a,b=ll=—s jV 

Similarly, the 0(mL) contribution to trC 2 can be calculated from the m = T = limit of 
the matrix C Q : 

.~ 2 . 1 ^^^sin^ 2-irla 2nlb 2nja 2-nkb . A . 

(&n)ik = "tttt > / r sin — — sm — — sin — —sin — — — . 24 

^ o>l« jy2 Z-j Z_/ gm *r& AT AT AT iV V ; 

a,o=l /=1 

The 0(mL) term of (C 2 )^ comes from zero modes pieces of C 2 where either a = in the 
{(f>4>) sum or a = in the (tttt) sum: 

o t , N-l s 2tt/6 „,,_ 27rfcb 

3mL coth (£M V V COS ^ COS ^ (25) 

0=1 (=— s jv 

(For (C e )jk, the indices have the range — s < j,k < s, while for (C )jk, we restrict to 



1 < i) & < s.) In the appendix, we describe how to perform the sums (23), (24), and (25) 
along with (|18[) in the the large N limit with s/iV held fixed. 



4 Raising the Temperature 

We present three arguments that the entanglement entropy depends exponentially on the 
ratio m/T in the limit T <C m. The first argument is heuristic and relies on the structure 
of the matrix C 2 . The second argument is based on our earlier calculation of trC 2 . The 
third argument is based on numerical evidence. We would like to show two things. The 
first is that for a fixed interval A, 

S(T) - S(0) ~ e~ m l T . (26) 

The second is that for two complementary intervals A and A = B, 

S A -S B ~ e- m ' T . (27) 

The first argument relies on the fact that the temperature dependence of C 2 comes 
entirely from the factors of coth(a; a /2T) in {(jxj)) and (tttt). The frequency uj a is bounded 
below by m. Thus we conclude that 

coth (g) < coth (g) = 1 + 1e~ mlT + 0(e-^' T ) . (28) 
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In other words, the matrix C has a low temperature expansion of the form 

C(T) = £7(0) + e- m/T 5C + ... (29) 

where the ellipsis denotes terms that are more exponentially suppressed. Now if C(T) 
has such an expansion, then the eigenvalues Pk(T) = ^(0) + e~ m l T bv\. + • • • will as well. 
Assuming Vk(0) -l/2> e _m//T , expanding eq. in the small T limit, one concludes that 
the entanglement entropy for a single interval shifts by an amount 



5S = 2 + 1/2) - HMV - l/2)}5vke~ m/T + . . . 



(30) 



implying the scaling (26). Assuming SS is different for an interval and its complement, one 



also concludes the scaling (27). 



While, the numerical evidence we present below suggests both scalings (26) and (27) 



are correct, there are some loop holes in our argument. An obvious problem is that the 
e -m/T t erm i n the small T expansion may vanish; the temperature dependence may be of 
the form e M l T for some M > m. A more subtle loop hole involves the fact that many of the 
ffc(0) are close to 1/2. In this case, the correction to the entanglement entropy SS can scale 
as (m/T)e -m / T instead of just e~ m / T . Numerically, we see no evidence for this behavior. 
Instead, in these cases we find that the logarithmic enhancement is not enough to make up 
for the smallness of Sv^\ these eigenvalues contribute negligibly to the entanglement entropy. 



The second argument for the scalings ( 26 ) and ( 27 ) is based on using tr C 2 as an estimate 



of the largest eigenvalue A e . Using trC 2 , we estimate the contribution of A e to S and infer 
the scalings from this contribution. The temperature dependence of a single interval comes 



principally from the leading coth(m/2T)/mL term in (14). One finds agreement with (26): 



S(T)-S(0)\x e ~e- m/T 



(31) 



Next we consider the entanglement difference Sb—Sa- This type of temperature dependence 
comes from the r 2 csch 2 (m/2T) piece of (14). One finds agreement with (27): 

7r mL 



S B (T) - S A (T)L . 



(1 - 2r) e~ m/T 



(32) 



VA,A e ,B 2 log TV 

We should emphasize that using trC 2 and the largest eigenvalue A e to estimate the 
temperature scalings is flawed. An obvious limitation is that we only have a result for tr C 2 
in the limit mL <C 1 while we expect the temperature scalings to hold more generally. 
A less obvious limitation is that despite the fact that A e is much larger than the other 
eigenvalues in the small mass limit, the logarithms in Q play a democratizing role and let 
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Figure 3: The contribution to 5S = Sa — Sb for the ten largest pairs of eigenvalues \j t A and 
Xj,B, arranged from largest to smallest. In this plot mL = 0.02, m/T = 1 and IjL = 1/5 
for region B. Note that the odd parity eigenvalues do not contribute. (N = 200 was used 
for this plot.) 



smaller eigenvalues contribute substantially to the entanglement entropy. For example, in 
this small mass limit numerical analysis shows that the dominant contribution to Sa — Sb 
comes from the second largest even eigenvalue (see figure [3]) . 



Our most convincing evidence for the scalings (26) and (27) is numerical and is presented 



in figures [4] and [5| Figure [4] demonstrates unambiguous evidence for ( 26 ) , not only for 
mL <C 1 but also for mL > 1. Figure [5^ displays unambiguous evidence for (27), again both 
for small and large values of mL. More ambitiously, we can try to investigate numerically 
whether the mL(l — 2r)/\ogN behavior of eq. (32) is correct as well. Figure [5^ provides 
evidence for the mL scaling. Figure [5|3 provides some limited evidence for the 1 — 2r behavior 
for large values of mL and for intervals with r ~ 1/2. However, we find no evidence for the 



logiV behavior of (32). 



5 Discussion 

As mentioned in the introduction, the original motivation for this paper came from the 
Ryu-Takayanagi proposal [TT] for computing the entanglement entropy of field theories with 
holographic dual classical gravity descriptions. In their proposal, the field theory lives on 
the boundary of the space-time in the dual description. Let C be the curve that separates 
region A from region B in the field theory. Let C also be the boundary of a minimal 
surface M that falls into the space-time. The proposal is that the entanglement entropy is 
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a) t b) t 

Figure 4: A log plot of the entanglement entropy 5S = S(T)—S(0) vs. m/T with an interval 
size £/L = 3/10. The points are numerically computed, and the line log(<55) = —m/T is a 
guide to the eye: a) mL = 5 x 1CP 3 ; b) mL = 5. (For both plots, the points were computed 
with N = 50, 100, 200, and 400. The data points for different values of N all lie roughly 
on top of each other.) 




Figure 5: a) A log plot of the entanglement entropy difference 5S = Sa — Sb vs. m/T for 
mL = 5 and 5 x 10 -3 , and an interval B of size 1/ L = 1/5. At fixed m/T, the larger mass 
points lie below the smaller ones. The line log(5S/mL) = —m/T is a guide to the eye. (The 
lattice was taken to have size = 200, but there is no noticeable difference between this 
graph and a graph with N = 100.) b) The entanglement entropy difference 5S vs. £/L for 
(from bottom to top) mL = 5 x 10~ 3 , 2, and 5. The mass to temperature ratio is m/T = 10. 
The line e m / T 5S/mL = 3m /T - 3/2 is a guide to the eye. (The lattice was taken to have 
N = 400, but there is no difference between this graph and a graph with N = 200.) 
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Figure 6: The two minimal surfaces Ma and Mb corresponding to a region A and its 
complement B when the dual space time contains a black hole (BH). 



proportional to the area of M: 



Area(M) 

SA = AGn ' (33) 



where Gn is Newton's constant. Assuming a unique such M, the entanglement entropy of 
a region and its complement are always equal, Sa = Sb- 

When the space-time contains a black hole, Ryu-Takayanagi modified their proposal to 
account for the existence of two minimal surfaces Ma and Mb- The entanglement entropy 
for A must be computed from the surface Ma that is deformable into A. Correspondingly, 
for region B, we must use Mb- For large black holes, Area(Mx) — Area(Ms) will come 
mostly from the differing amount of black hole horizon area that the two surfaces wrap (see 
figure [6]) . The Hawking temperature of the black hole corresponds to the temperature of 
the field theory, and thus this modification of the proposal provides a way for Sa — Sb to 
be nonzero for certain thermal field theories. 

However, there are instances where field theories at T > have dual gravity descriptions 
without a black hole. A classic example is the large N, strong coupling limit of maximally 
supersymmetric SU(N) Yang-Mills theory on S 3 x S 1 [12] . At temperatures small compared 
to the inverse radius of the 5 3 , the dual description is thermal AdS§ x S 5 . At a critical 
temperature T c , the gravity description undergoes a first order Hawking-Page phase tran- 
sition to a state with a large black hole. For the field theory, this transition is understood 
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as a deconfinement phase transition. 

On the one hand, their proposal implies that the entanglement entropy will serve as 
an order parameter for the phase transition: for T < T c , Sa — Sb = 0, while for T > T c , 
Sa — Sb / 0. On the other, at any finite N, we have a system at finite volume for 
which there can be no phase transitions. The transition from Sa = Sb at T = to 
Sa Sb at T > T c must be smooth. We conclude that the Ryu-Takayanagi formula is 
only valid in the strict large N limit, but it would be nice to understand the form of the 
1/N corrections. In principle, one should be able to compute the entanglement entropy for 
maximally supersymmetric Yang-Mills at weak coupling. In practice, such a computation 
is substantially more difficult, and we instead considered a 1+1 dimensional massive scalar 
field on a circle at T > 0. Morally, the regime T < m should correspond to the confining 
regime of the Yang-Mills theory where the fields get a mass through their coupling to 
the curvature of the S 3 . For our scalar field, we argued that in the regime T <^ m, the 
entanglement entropy difference scales as 

S B - S A ~ e~ m ' T • 

We conjecture that this type of scaling should be a generic feature of all gapped systems. 
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A Computing Traces in the Continuum Limit 

The basic technique used in computing trCg, trC^ and {tp\C^\ip) in the continuum limit 
N — > oo with n/N held fixed was to replace sums with integrals. However, there are three 
wrinkles in this procedure, two of which have already been hinted at in the text. The first 
is that we were not able to perform the integrals obtained by taking the continuum limit 
of the mode sums over a and b. Thus, we first performed the mode sums over a and b 
explicitly yielding sums over cotangents. For example, performing the mode sums for eq. 
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(34) 



while performing the mode sums for eq. ( 24 ) gives 

1 s 
AT2 E [' 



(C 2 ) 



i \-2 L cot + * " 1/2) " cot ]v~ (j + s + 1/2)+ 
i=l 



+ cot ^(j - I + 1/2) - cot £(j - s - 1/2) 
cot ^(k + I - 1/2) + cot ^(k - I + 1/2) 



(35) 



The second wrinkle is that naive integral approximations of the cotangent sums often 
include singular regions. Our strategy in this case was to add and subtract a sum that we 
could perform analytically but whose integral approximation had the same singular region. 



This procedure was already used in the text to perform the sum (21). The third wrinkle 
is that the integral approximations of the cotangent sums were often difficult to perform. 
Changing variables and using discrete symmetries reduced the integrals to known results in 
most cases. However, in two cases, we had to perform an integral we could not find in the 
books. 



Let us first sketch the computation of trC^, i.e. the trace of (35). Several of the terms 
in the sum have the structure 



1 s 

J ±± = Jp E cot jj( k ± - V2)) cot ^(k± (j - 1/2)) 



(36) 



k,j=l 



To perform these sums, we make the change of variables x = k + j and y = k — j. Using the 



same technique in eq. (21 ) to regularize the singular regions of the integral approximations, 
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one straightforwardly finds 

/++ = - 
J__ = s 



s z 1 



N 2 



, 2iVtan(7rs/iV) 
In v 7 ; +I + 7 



7T 



AT2 



7T 



In V ' ' +1+7 



7T 



+ 0(l/iV) , 
+ 0(l/iV 2 ) , 



/+- = I- 



1 



+ O (log N/N) . 



(37) 
(38) 
(39) 



An intermediate result necessary for the computation of I-\ is 



v E^ 1 )" cot £(* " V2) = " \ + *=*L cot ( ^ ) + 0(1/A^ 



AT 



k=l 



2N 



im 
N 



(40) 



The remaining pieces of tr C 2 can be rearranged in the following way 



2 E (cot^(& + i-l/2)+cot^(fc-i + l/2 



fc,i=i 



iV 



x (cot^(fc- s- 1/2) +cot-^(/c + s + l/2) 



2s \ 2s 

^cot^(y + l/2) -^ co t 2 ^(y + l/2) 



v 2/=0 
2s 2s 



y=0 



+ 2 E E cot ^(y + 1/2) cot ^(x + 1/2) . 



j/=0 x=2s—y 



N 



(41) 



The first sum on the r.h.s. of eq. (41) we performed in (21). The second sum can be 



performed using the same techniques: 

1 n 1 

^£cot 2 ^-l/2) = - + O(l/A0. 



k=l 



(42) 



The third sum requires more work and reduces to one of the two integrals we could not 
find in tables. Up to log N/N corrections, we may replace the third sum by the following 
integral: 

Hb) = / / cot(bx)cot(by)dydx , (43) 

JO Jl-x 

where 

A^E E cot^ (y + l/2)cot^(x + l/2)=(|) 2 /(^)+0(logiV/Ar). (44) 

y =0x=2s-y \ / \ / 

The integral over dy is trivial: 



1 f /i \ i sin(6) 

I ( b ) = v / cot (6a;) log . r^dx . 

o Jo sm(6(l - x)) 



(45) 
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We find that I'(b)b + 21(b) = -1 and that in the small b limit 1(b) = ir 2 /6b 2 + O(l). From 
these two facts, we deduce thatj^] 

m - £ - 1 . («) 

The quanties tr C 2 and (i/j\C 2 \ip) may be computed in an analogous way. As can be seen 



in eq. (34), there was one mode sum we were forced to do in the continuum limit: 

£ E = +fc ^) + op/*") ■ (47 > 

6=1 ' v 7 

All of the other mode sums we were able to perform explicitly. The remaining sums over 
cotangents are similar to cases treated above. We spare the reader almost all of the remain- 
ing details. In the computation of (ifj\C 2 \ifj) , we came across a second novel integral: 

! sin(6(l + x)) 12 



J(b) 



log- 



dx . (48) 



sin(6(l — x)) 

Similar to the strategy in computing 1(b), we find that J"(b) + 2J'(b)/b = —8 and that in 
the small b limit J(b) = n 2 /3 + 0(b 2 ). Thus we deduce that 

J(b)= l -(« 2 -Ab 2 ) . (49) 
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